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ABSTRACT ; 


Abstract  In  this  paper,  the  simulation  mode!  for  the  mutual  influence  between  the 
motion  of  the  rotational  parts  and  the  satellite  attitude  vibration  is  presented  by  using 
Lagrange  equation,  and  the  satellite  attitude  is  calculated.  The  influence  on  the  CCD 
image  geometric  quality  is  discussed. 

Subject  Term  Rigid  wheel  Satellite  attitude  Vibration  Mathematical  model 
Image  quality  Analyzing 


I .  Foreword 

Along  with  the  progress  of  economic  and  defense  programs  as 
well  as  scientific  undertakings,  requirements  of  geometric 
quality  of  satellite  remote  sensing  pictures  (such  as  resolution 
power,  in  regard  to  geometric  anomalies  and  positioning 
precision)  are  higher  and  higher.  Attitude  vibration  of  the 
satellite  is  one  of  important  factors  affecting  geometric 
quality,  especially  serious  in  the  case  of  remote  sensing  systems 
applying  linear  array  CCD  scan  imaging.  Therefore,  analysis  of 
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the  effect  of  satellite  attitude  on  picture  quality  is  one  of  the 
important  problems  concerning  designers. 

To  cope  with  the  continuous  progress  of  space  projects,  and 
the  more  and  more  complicated  dynamic  structure  of  spacecraft, 
the  craft  often  carry  rigid  rotating  components  and  elastic 
bodies  such  as  solar  sails.  This  article  analyzes  the  effect 
from  the  rotating  tape  in  a  recorder  on  the  geometric  quality 
(such  as  positioning  precision  and  internal  geometric  anomalies) 
of  imaging  with  respect  to  the  satellite-borne  CCD  multi-spectrum 
camera . 

Generally,  a  magnetic  tape  recorder  and  a  CCD  camera  are 
installed  at  the  bottom  (facing  the  Earth)  of  the  satellite.  The 
rotating  shaft  of  the  recorder  tape  is  parallel  to'  the  flight 
direction.  During  satellite  flight,  rotation  at  constant  linear 
speed  is  underway.  The  CCD  camera  applies  the  CCD  linear  array 
for  scan  imaging. 

During  installation  deviations  of  mass  center  and  shaft 
alignment  of  the  tape  recorder  wheel  may  occur.  In  the 
rotational  process  it  is  required  to  control  torque,  which 
controls  the  rotational  speed.  These  actions  may  cause  attitude 
vibration  of  the  satellite  and  directional  deviation  of  the  CCD 
camera,  thus  affecting  geometric  quality  of  CCD  imaging. 

In  this  article,  the  satellite  and  its  rotational  components 
are  considered  as  a  multiple  rigid  body.  By  applying  a  Lagrange 
equation  in  quasi-coordinates  and  another  Lagrange  equation  with 
Q  for  coordinates  in  generalized  sense,  a  simulation  model  of  the 
interdependent  effect  between  tape  wheel  rotation  and  satellite 
attitude  vibration  is  established.  In  addition,  attitude 
vibration  of  the  satellite  is  computed  by  using  a  numerical 
computation  method. 
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By  using  an  imaging  model  of  the  satellite-borne  linear 
array  CCD  camera  [5],  the  effect  of  attitude  vibration  on  picture 
quality  is  analyzed,  to  estimate  the  geometric  quality  of  remote 
sensing  pictures  by  CCD  cameras.  Besides  analyzing  the  effect  on 
picture  quality  due  to  attitude  vibration,  the  geometric  quality 
of  remote  sensing  pictures  with  the  CCD  camera  is  estimated.  The 
model  can  also  be  used  to  analyze  effects  on  picture  quality  due 
to  orbital  motion.  Earth  rotation.  Earth  curvature.  Earth 
projection  and  other  aspects.  In  the  article  the  situation  is 
taken  into  consideration  that  there  are  two  tape  wound  wheels  for 
each  tape  recorder. 

II.  Coordinate  System  and  Signal  Table 
The  Coordinate  System 

1.  The  coordinate  system  of  the  principal  inertia 

axis  of  the  recorder  wheel:  origin  0^  is  mass  center  of  the 
recorder  wheel;  O^X.^,  and  are  three  principal  inertia 

axes  of  the  wheel .  The  directions  of  the  three  axes  are  close  to 
the  three  respective  axes  of  OXiYjZi. 

2.  The  coordinate  system  OsXsYjZ,  of  the  principal  inertia 
axis  of  the  satellite:  origin  0^  is  the  mass  center  of  the 
satellite;  O^Xj,  O^Ys  and  OsZg,  respectively,  are  three  principal 
inertia  axes  of  the  satellite.  The  directions  of  the  three  axes 
are  close  to  the  three  respective  axes  of  OX2Y2Z2. 

3.  The  coordinate  system  OXiYjZi  of  the  recorder  wheel 
rotating  axis:  origin  0  is  the  combined  mass  center  of  the 
satellite;  OZi  is  parallel  to  the  shaft  axis  and  is  rigidly 
connected  to  the  recorder  wheel  (pointing  in  the  flight 
direction);  OYi  points  to  the  direction  of  Earth's  center  while 
OXj,  OYi  and  OZ^  form  a  right-handed  system. 
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4 .  The  coordinate  system  OX2Y2Z2  of  the  satellite  rotating 
axis;  origin  0  is  the  satellite  combined  mass  center;  OZ2  is 
parallel  to  the  shaft  axis  (pointing  in  the  flight  direction)  and 
rigidly  connected  to  the  satellite;  OYj  points  to  the  direction 
of  the  Earth's  center  while  OX2/  OY2  and  OZ2  form  a  right-handed 
system . 
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Fig.  1  Diagram  Showing  a  Model 
of  Magnetic  Tape  Recorder  Wheel 

Key;  1.  Main  shaft  of  recorder 
wheel;  2.  Bearing  axis; 

3 .  Foot  of  perpendicular 


Table  of  Symbols 

Superscript  and  subscript  r=a,  b,  c  or  d,  indicate  four 
recorder  wheels  respectively; 

- horizontal  direction  inertia  of  recorder  wheel; 

- rotating  shaft  inertia  of  recorder  wheel  (under  ideal 

situation  as  consistent  with  the  shaft  axis); 

^  -  dip  angle  of  recorder  wheel  A  principal  axis  relative  to 
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the  shaft  axis; 

.^,1  >0'^^  -  projection  of  3“  in  OXiYjZi; 

Jg  -  rotation  inertia  of  satellite; 


'  S 

'z^z 

A-, 

- projection  of  Jg  in  OX2Y2Z2; 


itir  - mass  of  recorder  wheel; 

mg  -  satellite  mass; 

^  -  vertical  vector  from  satellite  mass  center  to  shaft 

axis ; 

- projection  of  £  in  OX2Y2Z2; 

r  -  vector  from  recorder  wheel  A  mass  center  to  satellite 

mass  center; 

r°  -  vector  of  perpendicular  foot  on  shaft  axis  from 

recorder  wheel  A  mass  center  to  satellite  mass  center; 

- projection  of  r°  in  OXiY^Zj 

- control  torque  of  satellite; 

-  satellite  angular  velocity  in  inertia  space; 

( 6)352 / toy2 / <^z2 ) ^  projection  of  6)  in  OX2Y2Z2; 

Qr  -  relative  angular  velocity  of  various  recorder  wheels 

relative  to  satellite; 

(0,0,fir)^  - projection  of  S  in  OXiYiZj; 

T  -  total  kinetic  energy  of  satellite; 

$r=fir 

i(r - rolling  angle;  0 - elevation  angle;  and  (j>  - drift 

angle . 
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III.  Establishing  the  Model 

3 . 1  Model  of  the  recorder  wheel 

In  this  article,  in  the  satellite  are  installed  two  tape 
recorders  with  two  wheels  in  each  recorder,  a  total  of  four 
recorder  wheels,  marked  as  A,  B  and  C,  D  wheels.  A  recorder 
wheel  is  composed  of  a  disk  wound  with  tape.  In  the  photographic 
process,  the  tape  of  each  recorder  is  wound  around  another  disk 
(A-»B,  C-»D)  while  the  linear  velocity  Vq  is  constant. 

( 1 )  Main  rotation  inertia  of  the  recorder  wheel 

J'  +  J, 

Jp  is  rotation  inertia  of  the  recorder  wheel  while  Jq  is  rotation 
inertia  of  the  tape.  Jq  can  be  calculated  from  the  following 
equation: 

J.,  =  Y  +  Rl)Mr 
=  jMARl  -  Rl  +  /iV3) 


In  the  equation, 

Mr  is  mass  of  tape  wound  around  the  disk; 

Rp  is  internal  radius  of  the  recorder  disk; 

Rr  is  radius  after  tape  is  wound  around  the  disk;  and 
h  is  width  of  the  tape . 

( 2 )  Angular  velocity  and  radius  of  recorder  wheel 

When  starting  operation  all  tapes  are  wound  around  wheel  A 
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and  wheel  C.  With  simultaneous  start  and  reverse  rotating 
direction,  finally  all  tapes  are  wound  around  wheel  B  and  wheel  D 
as  the  following: 

a=-(2.,  <1>.  =-<!>.,  <P,=-<P, 

K  =  R^  =  R, 

In  brief,  it  is  assumed  that  for  each  loop  of  wheel 
rotation,  the  wheel  radius  adds  or  reduces  a  layer  of  tape 
thickness  H=0. 0295mm;  the  angular  velocity  and  angular 
acceleration  are  calculated  from  the  following  equation 


a  =  VJR^ 

=  [}\/(.Rr  ±  Hy  —  (subscript  r=a,  b,  c  or  d)  (1) 


(3)  Mass  of  tape  wound,  around  the  recorder  wheel 

=  M  ~  p  ’  V\  •  t 

Ml,  =  Mj  —  p  .  ,  I 

In  the  equations,  p  is  mass  of  a  unit  length  of  tape;  M  is  total 
tape  mass  wound  around  a  pair  of  recorder  wheels. 

3.2  Satellite  attitude  simulation  model 

Since  the  rotating  axis  of  satellite  rotation  components  is 
not  aligned  with  the  principal  inertia  axis  of  the  satellite,  and 
swing  is  caused  by  the  deviating  rotation  axis  from  mass  center 
position  of  the  rotating  components,  attitude  vibration  of  the 
satellite  will  be  induced. 

In  this  section,  a  Lagrange  equation  is  applied  to  derive  an 
attitude  model  of  the  three-axis  stabilized  satellite  with  four 
recorder  wheels.  Therefore,  the  equation  can  be  used  to 
calculate  the  effect  on  satellite  attitude  by  recorder  wheel 
rotation . 
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When  establishing  the  mathematical  model,  the  following 
aspects  are  assumed; 

1.  Both  satellite  and  recorder  wheel  are  rigid,  and  shafts 
are  rigidly  connected; 

2.  The  axis  of  rotatation  of  the  recorder  wheel  is  inclined 
while  installation  position  deviates. 


( 1 )  Kinetic  equation  of  attitude  motion 


By  using  Lagrange  equation  of  e  in  quasi-coordinates,  and 
another  Lagrange  equation  with  (j)  as  coordinates  in  the 
generalized  sense,  generalized  motion  equations  of  the  satellite 


as  follows; 

d  ■ 

IL ' 

—  <^r2 

jr 

dU.z  ! 

dui^Z 

d  ! 

1L\ 

- 

31' 

-j-  W.2 

<rr 

dt\ 

3u,y,  ! 

±i 

,71’  ■■ 

jr 

dt  1 

dUi,,  1 

-  ^.v2 

d(x)^2 

^yZ 

-  L,.  1 


(3) 


To  derive  the  kinetic  energy,  we  substitute  into  equation 
(3).  Following  linearization,  attitude  kinetic  equations  used 
for  simulation  analysis  are  derived. 


r  /* 

xZrZ 

-i_ 

1 

^X2X2 

-f 

‘^x2x2 

+ 

x2x2 

+  ^izxZ 

+  M(r%y 

y  3(r‘U  - 

e,2)']“x; 

2  + 

C^J-2v2 

+ 

J'lzyZ 

4- 

^  xZyZ 

+ 

J\zyZ 

+  ^izyZ 

+  3^^>z  — 

rU^irlz  - 

«v2  )  ]"v2 

4- 

C'^'x2j2 

+ 

+ 

^  xZzZ 

4- 

•^xZzZ 

-r  Jizxz 

4-  f^rlzi^xz 

= 

L.Z- 

C‘^x2x2 

-  Mrtzr% 

^^xZxZ^h 

JxZiZ^, 

3izzZ^J 

y  x2x2 

_L 

«^x2x2 

+ 

^"xZxZ 

4- 

^  x2x2 

]^^x2  — 

yx2^-2  +  ^UyZ 

^%ZyZ  +■  jizyz'l‘^yZ  ~  — 

^b')  *^x2i2^^=2  )  j xZzZ^^zZ  t  J  )  "i" 

f^xz(^:Z  i"  '^<i)^x2  ‘^/2^x2)(^y2  ®,2 ) 

//[ajx2(r°2  —  e.z)  +  “  ^?2)]^v2  + 

^x2C'^v2x2  f^’yZ^’izI^x  “f" 

+  JyZsZ^c  "i  J^yZsZ^t  '^.»2x2'^<^l 
“»,2[-^x2x2  +  f^r%(rlz  —  £,2)  +  MrlzUxz  —  rlz)^^2,  — 

^,2  L'^*2!2'^»  ”i  JzZsZ^c  “I”  J zZiZ^'i^ 
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iJ’yZyZ  +  J°yZyl  +  Aly'  +  '^kv2  +  ^  “  '■°2)']‘^.v2  + 

+  r,za  +  Auz  +  ^^2.,:  +  JU  +  ^'-'2(^.2  -  ''>2)  >.2  = 

L^,  —  [J”2.2  ~  IJ-r%rll']^a  —  JyZzZ^h  ~  ~  J^l^d  — 

[_j°xZy2  +  >j^*2.v2  +  'j^k*2  +  '^i2.»2l“'.r2  ~  {.^"ylyZ  +  ^ yZy^ 

^'yZyZ  >^v2.v2]“'.«2  ^  yZtl^'^tZ  +  '^a ) 

J*2,2(£t»,2  4-  iib)  —  ^‘ytzZ^^rl  +  ~  ^UzZ^^'-Z 

/^r"2(a>.2  +  Q.Vrlz  —  A2(ft2a2^^v2  —  “'.v2''x2)(£x2  —  2'x2  )  + 

—  £*2)  4"  ‘^yzC^xZ  ^i2^3''x: 

ai,2[Js2:2  ~  /^X2''«2)^a  +  JiztZ^b  +  J‘i2z2^c  +  J x2z2^<r\  + 

a.,2[^:2x2  4-  f^^zirlz  -  £,2)  -  t^lzi^zz  -  ^x2)Pa  4- 

<^x2V^2z2^ib  4-  J‘..2z2^z  +  ^^2-2^] 

[Ja2x2  +  J‘x2z2  +  J'xzzz  4-  J'x2z2  +  J^Z  +  -  r“2)]“'x2  4- 

[^:v2x2  +  J:.ZzZ  +  JyZ-Z  +  JyZzZ  +  JU  4-  f^lzi^yZ  "  ^°2)]‘^,2  + 

[^l2x2  +  J‘z2.2  4-  r.2zZ  +  J'xZzZ  +  J^zZ  ^  “  ''”2)^  4*  /'(Cx:  “  r^2)^]‘-.2  = 

L:,  -  [.J‘2-.2  -  -  r»2)  -  r»2(£x2  -  r°2)]-<^a  - 

J*2r2A,  -  J\2z~A  -  “  V^Zzl  4-  -^‘2.2  4"  J'x2=2  4-  ji2x2]‘^x2  “ 

[.>;2-2  +  ^"yZzl  +  ^'yZzZ  4-  j'‘y2z2']^yZ  "  ^^2x2  (“'-2  +  )  “ 

jx2-2(^x2  4-  a)  -  j‘y2z2i^z2  +  A)  “  ji2z2<^<^z2  4*  ^2,)  + 

t^i^yZ  -  ''v2)(^-2  +  -  £x2)(‘^x2  +  ^z')r\2  + 

fjt\_u>^2r\2  4-  (‘i-.zfyi  —  (<"=2  4-  •f2x)r“2 

IJ.[_(<ti:-  4-  ^^x)2'x2  ~  “'x:2'x2  ~  “'x2Ex2!]2'x2  4” 

a),2[(4“2x2  —  4-  JUzz^h  4-  J’xnz^c  4-  Jxi-.i^d]  — 

m,z'iU'‘y2..Z  4-  A^°2''^2)^x  4-  JyZzZ^i  4-  J'yZzZ^c  4-  JUzZ^A 
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(2)  Satellite  kinetic  equation 


To  analyze  the  effect  on  imaging  by  attitude  motion,  it  is 
also  required  to  calculate  the  attitude  angle  of  the  satellite 
based  on  angular  velocity  e  of  the  satellite  relative  to  inertia 
space . 

Assume  that  Q*  is  the  orbital  angular  velocity  of  the 
satellite,  by  taking  small  angle  approximations  of  t,  4)  and  fl*, 
while  neglecting  high-order  small  quantities,  equations  of 
satellite  attitude  motion  are  obtained. 


99=  Wj.  —  O*  *  (p 

Wy  -  n* 

(pz=.  fl*  •  <p 


t 


(7) 


IV.  Simulation  Analysis 

4 . 1  Calculation  of  attitude  variation 

The  following  parameters  are  picked  out  in  simulation: 
Principal  rotation  inertia  of  the  tape  disk: 

=  0.  008kg.  .  j;  =  0.  008kg.  =  0.  01 62kg. 

r’ =  =  (2X10''m.  2X10'‘m,  0.  2m)'^ 

=  (0.  0167rad,  0.  0167rad.  0)'^ 
e=  (e^j,  0)“^  =  (0.  4tn,  —0.5m,  0)"^ 

Take  Qq  as  the  initial  velocity  of  recorder  wheel. 

By  using  equations  (4)  through  (7),  we  select  parameters  and 
apply  the  numerical  calculation  method  to  derive  variation  of 
satellite  attitude  (see  Fig.  2).  Then  from 

Li  —  +  </j,,,d£u,  /dt 

A,  =  U'liyzdWi  4-  J'ylizdui,  +  j;„,jda;J/di 
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we  derive  the  interference  torque.  See  Fig.  3  for  calculation 
results . 


Fig.  2.  Satellite  attitude  variation 
without  considering  start  process 


Fig.  3.  Interference  torque  due  to  rotation  of  magnetic 
tape  wheel  when  starting  process  is  not  considered 

Key:  1.  Torque;  2.  Torque  along  direction  of  Y  axis; 

3.  Torque  along  flight  direction 
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4.2  Imaging  quality  of  CCD  camera 


To  further  analyze  the  effect  on  CCD  imaging  quality  by  tape 
recorder  operation,  the  physical  process  of  imaging  is  applied  to 
establish  an  imaging  model.  By  using  this  model  and  the  above- 
mentioned  data  of  attitude  angle,  a  simulation  graph  is  plotted. 
Then  another  simulation  graph  is  plotted  by  using  attitude  angle 
at  zero  (in  the  case  without  recorder  wheel  rotation)  as  the 
standard  graph.  The  picture  amplitude  is  equivalent  to  5400  x 
5400  picture  elements.  A  point  is  picked  for  every  600  picture 
elements  in  each  line  and  each  row  to  be  used  as  the  control 
point  for  precision  check. 

( 1 )  Calculate  mean  square  deviation  of  the  corresponding  points 
for  two  pictures 

DDD=  -^0.)^  +  ^ 

In  the  equation,  (Xi,  Y^)  is  the  point  in  the  simulation  graph 
with  tape  wheel  rotation;  (Xoi,  Yoi)  is  the  point  in  the 
simulation  graph  without  tape  wheel  rotation. 

These  data  indicate  error  in  absolute  positioning  as 
affected  by  attitude.  See  Table  1  for  calculation  results. 

Table  1 .  CCD  imaging  indexes  of  simulation  calculation 

Condition  DDD  DDTM 

Considering  start  process  0.00608m  0.00376ra 

No  start  process  0.01219m  0.00686m 


(2)  Calculate  error  from  internal  anomaly 
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Mean  value  of  error  AX  =  —  V  (X,  —  X,.) 

Ay=lV(y,_y^), 

Standard  deviation  for  single-direction  error 

sx’=  ->^0,)'  -  (AX)*]/(«-  1) 

sy  =  -^0.)*  -  (Ay)*]/(«  -  1) 

standard  deviation  of  error 

DDTM  =  [(fiX)*  +  (Sy)']*'* 

DDTM  indicates  the  relative  positioning  error  due  to 
parameter  deviation,  thus  revealing  the  internal  anomaly  of  the 
graph.  See  Table  1  for  calculation  results. 

Consideration  is  given  to  the  starting  process;  that  is, 
when  calculating  attitude,  the  initial  velocity  of  the  recorder 
wheel  is  0  with  consideration  of  the  start  process. 

When  t  <  6s, 

=-  da  =  V„/(r^  X  6) 
da  =-  da  =  Vo/(r^  X  6) 

Without  starting  process,  that  is,  during  calculations  of 
the  attitude,  the  initial  velocity  of  tape  wheel  is 
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